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abstract 

While the literature abounds with sources for 
either near critical or high pressure heat transfer 
data, aside from fluid hydrogen, a dearth of informa- 
tion exists on a unified experiment over a wide range 
of reduced temperature and reduced pressure. Forced 
convection heat transfer data for fluid nitrogen 
flowing through a vertical tube were acquired for 
0.87 < Tj. < 1.5 and 1.2 < < 5, (P^ = 3.A2 MPa, 

Tj. = 126,3 K). For design purposes, these data were 
assumed to simulate the heat transfer characteristics 
of LOX (P^ * 5.08 MPa, T^, = 154.8 K) . Comparisons 
between LN 2 and LOX data over the pressure range for 
Ty < 1 using a modified form of Petukhov's correlation 
factor wad sho.n, qualitatively, to substantiate the 
simulation. Other parameters (buoyancy, Stanton mrni- 
ber, and relative velocity) , were found to group 
those data; assessment of these parametric inter- 
relations must be made before the simulation can be 
termed quantitative. 

NOMENCLATURE 

Cp specific heat, J/g-K 

0 diameter, cm 

F thermophysical properties parameter, eq. (2) 

Fr Froude number 

fj^-f^ function forms - see text for definition 
G mass flux, g/cm^-sec 

g body force, N 

H enthalpy, J/g 

h heiVt transfer coefficient, W/cm^-K 

k thermal conductivity, W/cm-K 

L reference length, cm 

Nu Nusselt number 

NUj. reduced Nusselt number, Nu^j^/Nu^^^c 

P pressure 

Pr Prandtl number 

q heat flux, W/cm^ 

R gas constant, (MPa) (cm^)/ (g) (K) 

Ra Rayleigh number 

Re Reynolds number, Gd/p 

Sr Sterroann Parameter 

STAR Category 34 


T temperature, K 

T* transposed critical temperature, K 

V Specific volume, cm^/g 

W mass flux, g/sec 

X axial position, cm 

a thermal diffusivity, cm^/sec 

8 volume expansitlvity, (3 In V/3T)p, 

X latent heat of vaporization, J/g 

6 Dittus-Boelter properties parameter, 

p density, g/cra^ 

U viscosity, g/cm-sec 

Subscripts 
B bulk 

c critical 

calc calculated , 
exp experimental 

r reduced parameter (eg. Tj. = T/T^,) 

w wall 

1 heavy component 

2 light component 
Superscript 

* wall to bulk property ratio 

INTRODUCTION 

One of the proposals for a second generation 
shuttle vehicle suggests that the engine regenerative 
coolant be the oxidizer rather than a fuel as is usu- 
ally the case (1) . The feasibility of such a proposal 
depends heavily on the ability of oxygen to function 
as a coolant. While there exists a wealth of heating 
and cooling data in the literature, there exists a 
dearth of information on the heat transfer character- 
istics of liquid oxygen (LOX) in the near critical 
region at elevated pressures (proposed engine oper- 
ating pressures, 20 to 50 MPa). In this paper we pro- 
pose to make use of data from similar fluids, thus , 
tapping this source of information and obviating mahy 
of the safety problems associated with a heat transfer 
test of fluid oxygen at high pressure. 

Recently, Hendricks and Simoneau(^) proposed and 
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demonstrated that two phase choked flow of fluid oxy- 
gen could be accurately simulated using fluid nitrogen. 
The principle of corresponding states was successfully 
applied to two phase choked flows of simple fluids 
over a broad range in reduced temperature and pressure 
as well as for a variety of nozzles, orifices, and 
slits. 

This same concept was assumed to apply to the 
heat transfer characteristics of liquid oxygen. As a 
result, two experimental heated tube programs were 
initiated, one at Lewis Research Center using fluid 
nitrogen to simulate LOX heat transfer and the other 
at Aerojet General using LOX to simulate actual engine 
operating conditions. 

The initial phase of the latter study was re- 
ported in reference 3. Rousar and Miller (3) pre- 
sented an equation for the cooling characteristics of 
fluid oxygen based on the modified property ratio pa- 
rameters correlation of Petukhov (4) for pressures to 
34.51 MPa over a range of experimental ccnditions 
simulating engine coolant requirements . 

The purpose of this paper is to consider simi- 
larities and data to support the thesis that fluid 
nitrogen simulation of fluid oxygen heat transfer 
characteristics is proper and can be justified. 

LIQUID NITROGEN APPARATUS AND INSTRUMENTATION 

The experimental test facility and apparatus is 
essentially that reported in reference 7 except that 
the working fluid is nitrogen. The joule heated test 
section is a 0.635 cm o.d. by 0.050 cm wall Inconel 
tube 43 cm long (L/D = 81) . The surface temperatures 
were measured and recorded at 19 stations along the 
tube. Three pressure taps in addition to the mixinc 
chamber conditions provided for pressure measurements. 
The accuracy of the recording system and transducers 
is discussed in reference 7. The primary measurements 
are heat flux, mass flow rate, surface temperature, 
and (inlet/outlet) conditions (i.e., the standard 
q-W-AT type experiment) . B'or most jf the data here- 
in, the system heat balances are within ±5 percent 
and redundant flow measurements compare within ±2 per- 
cent, 
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heated tube experiment with liquid oxygen in the near 
critical region. In 1961 Petukhov (4) presented a set 
of empirical property ratio parameters which appeared 
useful in correlating near critical heat transfer data 
These and many other techniques both theoretical and 
empirical for evaluating the heat and momentum trans- 
fer characteristics of near critical fluids are pre- 
sented in a survey article by Hendricks, Simoneau, and 
Smith (6) , 

From the literature survey (.§) it is quite appar- 
ent that we neither understand the fundamentals of 
heat and mass transfer in the near critical region uor 
can we recommend a universally empirical correlation 
for existing forced convection data. However, we can 
readily summarize some of the qualitative concepts re- 
lating to the heat transfer in the near critical re- 
gion (^) . For example, the modified Petukhov parame- 
ter (M and the Dittus-Boelter property parameter = 
c0.4j^0.6/|j0.4 evaluated at some reference temperature 
are useful in grouping near critical data. {Roth v;ill 
be used herein.) 

We know there exist rather large (L/D) effects 
Implying thermodynamic nonequilibrium and in some 
cases large volp-metrlc expansion effects which imply 
local and/or global instabilities. We also know tiiat 
near critical p-hydrogen data nan be correlated by a 
modified Dittus-Boelter type equation over a narrow 
range of reduced temperatures (0.8 to 1.5) at reduced 
pressures of 10 to 13. Furthermore, the fluid behav- 
ior of liquid and gaslike flows on either side of the 
critical region can be correlated using a Dittus- 
Boelter equation modified by some form of a (T^/T^) 
correction (^) . 

With these concepts in mind, we elected to exam- 
ine several correlating parameters which, afforded some 
degree of success in grouping near critical heat 
transfer data. One of these parameters was of the 
modified Petukhov type (,4,^), and quite similar to 
that presented by Rousar and Miller (3). With some 
adjustments of the exponents, this correlation group 
could Psirve as a common technique to examine the 
the'-xs that the heat transfer characteristics of fluid 
oxygen can be simulated using fluid nitrogen; this is 
more important than offering still another correlation 
technique, although some alternate parameters ate dis- 
cussed in appendix A. 


While the concept was to simulate fluid oxygen 
heat transfer at pressure levels characteristic of 
engine operating conditions, fluid nitrogen opera- 
tions were limited by the facility to 17.25 MPa. At 
Che same reduced pressure, this is equivalent to 
fluid oxygen operation at 26 MPa, or about half of 
the upper bound but well within the proposed oper- 
ating pressures. Some thermodynamic critical param- 
eters for fluid oxygen and nitrogen are given in 
table 1. 

The range of experimental parameters are 

Heat flux to 750 W/cm^ 

Pressure to 17 MPa (1.2 < < 5) 

Mass flow rate to 840 g/sec 

Wall temperature to 650 K 

Bulk temperature 110 to 200 K (0.87 < Tn/T. 

< i-6) 

DISOUSSiON AND RESULTS . , 

In 1957, Powell (^ reported the results of a 


Electing to leave the equation in a grouped for- 
mat to illustrate the significance of each, the form 
chosen for comparison becomes 
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For the fluid nitrogen data the value of Y 
grouped primarily between 0.65 and 1.15 with a possi- 
ble trend to lower values at low reduced enthalpy and 
higher at the higher reduced enthalpy. The fluid 
oxygen data of Powell (5) showed similar trends but 
group between 0.9 to 1.4 with a little more scatter. 
As the grouping of Y was similar to that of fig- 
ures 1 and 2 , one can deduce the data can be grouped 
in a similar way without altering the Reynolds nucber 
exponent . 


The first parameter set of equation (1) essentially 
Includes all of the variables in the heated tube ex- 
periment (experimental control group) . The second 
parameter set represents thermophysical property 
variations in response to the system operating condi- 
tions. The therraophyslcal property data were calr.u- 
lated using Subroutine GASP (7) . 

While we have investigated to a limited extent, 
the influence of the functional forms f^, f 2 » £ 3 , 
and f^, the most significant is the group. 


From figures 1 and 2' and the above discussion, 
one may conclude that the aodiEled Petukhov parameter 
(eq. ( 2 )) is effective in qualitatively grouping the 
experimental heat trans'er data for fluids oxygen and 
nitrogen in a similar r anner as proposed; however, 
there is considerable scatter in these results. 

Quantitatively, however, the oxygen and nitrogen 
Nusselt numbers are at different levels. The Nusselt 
number levels appear to be related through their T* 
ratios 


In order to compare the fluid nitrogen data with 
the fluid oxygen data of Rousar and Miller ( 6 ) which 
covered a large range of heat transfer conditions ns 
well as the data of Powell (^) , the functional forms 
f^ - f^ have to be modified as follows: 




where T* is given in appendix A, and p « 2 for 
equation (3) and p s 1 for equation (4). 

OTHER PARAMETERS 

Other parameters related to a rapid change in 
fluid properties that can influence heat transfer 
such as Gra/Jhof number, penetration length (Stanton 
number) and relative velocity effects were examined 
for similarities when applied to fluids oxygen and 
nitrogen. Discussion of these parameters is pre- 
sented in appendix A. In presenting the figures of 
this section it is useful to use reduced Nusselt num- 
ber, equations (la) and ( 2 ). 


Under these conditions, the fluid oxygen data of 
Powell (, 5 ) are presented in figure 1 as a function of 
reduced enthalpy. The Value of Y for these data 
group between 0.8 and l.i3. An inspection of the data 
presented by Rousar and, Miller indicates that the 
value of Y for most cf their data fall between 0.7 
and 1.3. The ±30 percent scatter bands on these and 
subsequent figures aid in comparison of the nitrogen 
and oxygen data sets . Considering the large range of 
parameters covered by the experiment of reference 3, 
the agreement with the Powell data is quite good. 
Using these same parameters, the fluid nitrogen data 
are plotted in figure 2. The value of Y for these 
data is between 0.4 and 0.8, which is qualitatively 
quite similar to that of oxygen. 

In order to determine if the increase in 
Reynolds number exponent from 0.8 to 0.95 (combining 
£4 of equation (3) with equation ( 2 )) sUppre.sses 
other effects such as (L/D) we assigned the to 

f^ functions as: 



To determine if the influence of body forces _ 
(volumetric expansion effects) is the same for 
fluids nitrogen and oxygen, the reduced Nusselt num- 
ber was plotted as a fufiction of the parameter 
8 AT/a^. Comparing figures 3(a) and (b) indicates a 
decrease in heat transfer as the parameter (Grashof 
effects) increases for both fluid oxygen and nitrogen. 
There also appears to be a threshold value below 
which the parameter has little effect on Nu^ or the 
liquid data simply group at large values of the pa- 
rameter. 

To determine if the Stanton number (penetration 
length) Influences heat transfer to fluid oxygen and 
nitrogen in a similar manner, the reduced Nusselt 
number was plotted as a function of (q/G(H^j - Hg)). 

The resulting plot showed trends which were similar 
for fluids oxygen and nitrogen and, furthermore, 
similar to the hydrogen data of reference 8 . These 
data were found to be nearly dependent on the square 
root of the reduced Nusselt number and to have a 
greater dependence on Prandtl number. 

Nu_ - - St3/2 or (5) 

«“calc Pr0*7 

tiding equation (5) , the oxygen data of Powell 
and nitrogen data are plotted as figures 4(a) and (b) , 
respectively. At low values of reduced enthalpy, the 
trends are similar, while at elevated values no com- 
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parison can be made although the oxygen data do appear 
to have a wavelike form. 

The final parameter investigated (1 + (3 In V/ 
3H)6H) indicates that reduced Nusselt number decreases 
as the relative (slip) velocity increases, see figures 
5(a) and (b) for fluids- oxygen and nitrogen, respec- 
tively, Such behavior is Indicative of an increased 
boundary layer thickness similar to that associated 
with film cooling. 

For each of the above parameters tested, fluid 
nitrogen and fluid ojQ^gen behave in a somewhat similar 
manner; however, there are several issues to be re- 
solved before the heat transfer behavior can be pre- 
dicted quantitatively, 

CONCLUSIONS 

Figures 1 and 2, and the agreement with the 
large data set of reference 3, demonstrate that ex- 
perimental heat transfer results for fluid nitrogen 
can be used qualitatively to simulate heat transfer 
to liquid oxygen. Quantitative agreement appears 
satisfactory provided the proper level or adjusting 
constant can be determined. This constant may be re- 
lated to the ratio of transposed critical tempera- 
tures . 

Figures 3 to 5 further demonstrate the qualita- 
tive agreement in the heat transfer behavior &f 
fluids oxygen an-d nitrogen; however, the parameters 
tested indicate dependence on terms such as buoyancy, 
Stanton number and relative velocity which are not 
part of the correlation used to produce figures 1 
and 2. Answers to such interrelations must be deter- 
mined before the simulation can be termed quantitative. 

It is suggested that heat transfer tests be con- 
ducted on fluids nitrogen and oxygen in the same fa- 
cility with the same test apparatus with test parame- 
ters spanning the same range of reduced conditions. 

Such a procedure which was so successful in the work 
of reference 2, should also prove productive here. 

APPENDIX A 

Near Critical Convection 

Many authors have proposed and found that dete- 
rioration of heat transfer in the critical region is 
body force dependevi. However, in the case of hydro- 
gen, the temperature profiles remained unaltered for 
vertical up or down flow in a 0.851 cm (0.335 In.) 
diameter tube at Reynolds numbers ranging from 1 to 
5 M. At these Reynolds numbers, the inertial forces 
are quite large and deterioration appeared to have an 
analogy to the transition from nucleate to film boil- 
ing. The nature of the axial temperature profile ap- 
pears to be a significant factor, in determining 
whether or not body forces are Important. In any 
case, the relative magnitude of the forces involved 
should be considered as the hydrogen data of refer- 
ence 7 are clearly separated into two bands. In the 
first band, are data from two larger diameter Lubes , 
1.113 and 1.28 cm (0.438 and 0.507 in.) and in the 
second band are the data of three smaller tubes 
0.478, 0.559, and 0.851 cm (0.188, 0.22, and 0.335 
in.) diameter. A Frdude number effect was used to 
combine these bands. Boussinesq (9) pointed out that 
for highly turbulent flows in natural convection, the 
Nusselt number should be related to the product 


Rayleigh x Prandtl: 


Nu = RaPr = p AT 
^2 


Here, the significant parameter is 8 AT/a as the 
test section of references 3, 5, and herein are nearly 
the same diameter with vertical upward flow. 

8 AT represents the body force (or volinne expan- 
sion) potential and a is the molecular diffusivity 
which is related to the periodic behavior of energy 
transport in the vicinity of the wall. In the near 
critical region, 8 is large, o is small, and in our 
case AT is large; consequently the parameter can be 
quite large. 

It should also be pointed out that the Froude num- 
ber relates to the body force potential, relative or 
slip velocity and penetration length when using the 
Boussinesq formalation (eq. (Al)), 


Fr ^ 


St'^^Nu 

RaPr 


8 AT = hr P 4X 


Stanton Number or Dimensionless Penetration Length 

The interchange of light and heavy components in 
the vicinity of the wall is analogous to a jet of fluid 
moving toward and away from the wall (^) . A ratio of 
the normal velocity to the free stream velocity can be 
considered to be a dimensionless penetration distance. 
The Stanton number can be interpreted to be a dimen- 
sionless penetration parameter. Stermann ( 10) defined 
a pseudo-Stanton number for two-phase fl; w as 


XpU ' 


. St = ^ 


Relative Velocity 

A form for the relative velocity , which may be 
viewed as a slip condition, can be postulated using 
the concept that the mass influx equals mass outflux 
of a control volume in the boundary layer or 

2 = influx ^ ^2^A (A4) 

^outflux 


A 1 /A 2 = -S given in referen =. 8, and 

V 2 /V 2 may be deter. Ined from a truncated; series re- 
lating volume and energy across the boundary layer as 


- 1 = 14 - 1 - 


and the relative velocity becomes 


A large value of this ratio would imply small volume 
toward the wall and large volume away from the wall. 

Transposed Critical Temperature 

While not precisely correct, the transposed criti- 
cal temperature for fluids nitrogen and oxygen may be 
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represented by 

+ 2s(|^y+ (V -23.5) (A7) 

.^'or 1 < (Pg/Pj.) < 4 And constant at (Pg/P ) = 4 

for higher pressures to perhaps 6Pc {^) . 
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TABLE I. - CRITICAL CONSTANTS FOR FLUID 
OXYGEN AND NITROGEN 


Fluid 

Pressure, 

MPa 

temper- 

ature, 

K 

Density, 

g/cm3 

Compres- 

sibility 

Nitrogen 

3.417 

126.3 

0.3105 

0.294 

Oxygen 

5.082 

154.8 

.4325 

, 292 
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Fig- 1 Influence of the functional group of equation 13) on the 
fluid oxygen heat transfer data of Powell (ref. 5) as a 
function of reduced enthalpy. 
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Fig. 2 Influence of the functional group of equation (3) on the 
fluid nitrogen data as a function of reduced enthalpy. 













